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I Abstract 

We introduce conformal coupling of the Standard Model Higgs field to gravity and 
discuss the subsequent modification of i? 2 -inflation. The main observation is a lower 
temperature of reheating which happens mostly through scalaron decays into gluons 
due to the conformal (trace) anomaly. This modifies all predictions of the original R 2 - 
q inflation. To the next-to-leading order in slow roll parameters we calculate amplitudes 

£3 and indices of scalar and tensor perturbations produced at inflation. The results are 

compared to the next-to-leading order predictions of i? -inflation with minimally cou- 
pled Higgs field and of Higgs-inflation. We discuss additional features in gravity wave 
J> signal that may help to distinguish the proposed variant of i? -inflation. Remarkably, 

the features are expected in the region available for study at future experiments like 
BBO and DECIGO. Finally, we check that (meta)stability of electroweak vacuum in 
the cosmological model is consistent with recent results of searches for the Higgs boson 



(N 



at LHC. 



1 Introduction and Summary 

The Starobinsky model of inflation [1] is the first, yet realistic example of new physics capable 
of solving major problems of the Hot Big Bang theory, see e.g. [2]. It exploits dynamics of 
gravity sector, modified by a quadratic in scalar curvature term added to the gravity action. 
The attractive feature of the model is that one and the same force-gravity-is responsible for 
both inflation and subsequent reheating of the early Universe. Such minimality is of some 
interest, given the absence of any direct evidence of relevant new physics in laboratory and 
accelerator experiments. 
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In this paper we consider the Starobinsky inflation with matter sector described by the 
Standard Model of particle physics (SM) which scalar sector is slightly modified. Namely, 
we add a conformal coupling of the SM Higgs field to gravity. This term leaves intact the low 
energy phenomenology of the SM, but impacts on the history of the early Universe. Indeed, 
we found that with the Higgs boson becoming conformal at high energies, reheating of the 
Universe takes place later and occurs via gluon production due to the conformal (trace) 
anomaly. The idea of conformal anomaly being responsible for reheating was discussed in 
literature, e.g., [3, 4]. Here it is natural consequence of the conformal symmetry in our 
model. 

Lower reheating temperature implies longer matter dominated stage between inflation 
and reheating. This modifies all predictions for power spectra of scalar and tensor pertur- 
bations generated at inflation. Likewise, this modifies predictions for gravity wave signals 
expected from nonlinear structure dynamics at post-inflationary stage. These are special 
signals in gravity waves given the long-lasting post-inflationary matter dominated stage. 
Remarkably, the signals fall in the region expected to be reached by proposed future experi- 
ments like BBO [5] and DECIGO [6] on searches for gravity waves. These signals have been 
proposed [7] as signatures of i? 2 -inflation. The same is true for our variant with conformal 
Higgs, where the features in gravity wave spectrum are expected at different frequencies, 
which allows to test the model. Finally, the non-minimal coupling provides with additional 
term in the Higgs effective potential, which becomes important at large scalar curvature. 
This can change the answer to the question: in which vacuum does the Higgs field fall in the 
expanding Universe, given the value of the Higgs self-coupling (or the Higgs boson mass). 

We address all these issues below. The model is presented in Sec. 2, and reheating 
is studied in Sec. 3. Predictions for amplitudes and spectral indices of scalar and tensor 
perturbations are obtained in Sec. 4. The calculations are performed to the next-to-leading 
order in slow roll parameters and the results are compared to similar predictions obtained 
there for the original Starobinsky model and for Higgs- inflation [8]. All the three models 
exhibit the same inflationary dynamics, so the only difference is in reheating temperature, 
which helps to distinguish the model predictions. In Sec. 5 we discuss the gravity wave 
signals expected in the model: one comes from inflation, others from nonlinear evolution of 
inhomogeneities at post-inflationary matter-dominated stage. Sec. 6 is devoted to analysis 
of stability of electroweak vacuum. There we estimate the lower bound on the Higgs boson 
mass, corresponding to the viable cosmological evolution in the model and find it to be 
consistent with recent results of LHC. 



2 



2 The model description 



The Starobinsky model of inflation is described in the Jordan frame by the following action 
[1, 9], 

S = ~^Y~ j d 4 X - ^—^j + Smatter ■ (1) 

Here the reduced Planck mass is M P = M P i/\^8tt = 2.4 x 10 18 GeV and S ma tter in our 
case refers to the SM action. At large values of scalar curvature R model (1) allows for 
an inflationary stage in a slow roll regime: a scalar degree of freedom (dubbed scalaron), 
emerging when i? 2 -term is added to the gravity action, plays inflaton. Scalaron quantum 
fluctuations evolving at inflationary stage freeze out with amplitude ~ /i when exit the 
horizon. Later they give rise to matter perturbations, which amplitude is fixed by the global 
fit to cosmological data, consequently the value of parameter \i equals [9] 

H = 1.3 x lCT 5 M P = 3.1 x 10 13 GeV . (2) 

After inflation the Universe expansion is driven by oscillating massive scalaron field re- 
sponsible for the effective matter-dominated post-inflationary stage. Later, scalaron oscil- 
lations decay into the SM Higgs bosons due to gravity interactions and a hot stage in the 
Universe starts with temperature [10] 

T r S = 3.1 x 10 9 GeV . (3) 

In this paper we introduce in model (1) conformal coupling of the SM Higgs field to 
gravity. Then the part of action with the Higgs doublet % reads (we omit irrelevant for our 
study Yukawa terms): 

S H = J^d A x(^RH ] H + D^D^H-^{H ] H-v 2 ) 2 ^j . (4) 

After the conformal (Weyl) transformation to the Einstein frame 

9r eVW*'"^ (5) 

action (1) takes form [11] 

r f m 2 1 \ 

S = J^-gd A x 1—fR + - V{4>)\ + S, 



matter j 



Here S matter is the conformally transformed action of matter (SM fields). Any conformal non- 
invariance in the matter sector produces interaction between scalaron <fi and SM particles. 
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3 Reheating via the conformal (gauge) anomaly 



In our model (1), (4) the SM Higgs is conformal at large values of the field and hence it 
decouples from scalaron. Then the strongest relevant coupling between scalaron and SM 
fields come from the conformal (gauge) anomaly which takes care of reheating in the model. 
Indeed, the conformal transformation (5) yields scalaron coupling to the trace of energy- 
momentum tensor of matter fields Tff: 

S lnt = Jv^d'x^TZ. (6) 

The relevant terms in (6) are due to conformal (gauge) anomaly, see e.g. [12]: 

Here F£ v stand for gauge field tensors and b a are the first coefficients of /3-functions for 
corresponding gauge coupling constants a; for the SM couplings of U(l)y, SU(2)\y and 
<S77(3) C gauge interactions these coefficients are y, — y and —7, respectively. 
As a result, the scalaron decay rate into the gauge bosons is 

_ b 2 a a 2 N &d] 

1^2 bosons- {.*) 

where A^j = 1, 3, 8 for U(l)y, SU(2)w and SU(3) C gauge interactions, correspondingly. 
Values of a must be taken at the scale of and we obtain them by making use of the 
numerical code [13] operating with the SM 3-loop /3-functions [14]. 

Scalaron decays mostly into gluons, which immediately rescatter producing all the SM 
particles 1 . The scalaron total decay rate in our model with conformal Higgs field (1), (4) 
is about 140 times lower than that in the model with the Higgs field minimally coupled to 
gravity (1). We define the reheating temperature of the Universe after inflation as temper- 
ature at the moment of equality between the energy densities of scalaron condensate and 
relativistic matter [10]. Then numerically 

T rch = 1.1 x (?; 1/4 (T rch ) = 1.4 x 10 s GeV , (9) 

where effective number of degrees of freedom in the plasma of SM particles is g* (T ren ) = 
106.75. 

1 Amplitudes of scalaron direct decays to pair of SM massive fermions are suppressed by corresponding 
Yukawa coupling constants and ratio of the Higgs field ~ H, see Sec. 6, to the Planck scale. Rates of three 
and four body scalaron decays to SM particles are strongly suppressed by coupling constants and the phase 
space volume. 
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4 Parameters of scalar and tensor perturbations pro- 
duced at the inflationary stage 

In both variants of i? 2 -inflation (with ordinary and with conformal Higgs field) the scalaron 
potential at inflation is the same. Therefore, the only difference in predictions for parameters 
of scalar and tensor perturbations generated at inflation is due to different number of e- 
foldings because of different reheating temperature in the models (for detailed explanation 
see e.g. [7]). Since parameters of perturbation power spectra depend on the reheating 
temperature very mildly (logarithmically) and the latter differs in the two models only by 
factor of 20, cf. Eqs. (3) and (9), they must be evaluated to the next-to-leading order in 
slow-roll parameters. 

The procedure is described e.g. in Refs. [15, 16]. For the leading order estimates one usu- 
ally exploits the number iV e of e-foldings passed after the perturbation of a given conformal 
moment k exited the horizon [2, 17], 

Ar , fa(k)\ n 1, /1.4xl0 8 GeV\ , , 

N e = log l±±\ « 53.27 - - log ( — J . (10) 

Here a(k) and a e refer to the scale factor at the moment of horizon exit and at the end of 
inflation, respectively. The value of k is chosen to match the WMAP pivot scale k/a = 
0.002 Mpc" 1 , where ao is the present scale factor. 

For the next-to-leading order calculations a more convenient measure of the moment of 
horizon exit is [16] 

N, - log [«) , (11) 



where H{k) and H e stand for the Hubble parameter at the exit and at the end of inflation. 
The latter is defined as the moment when the Universe stops to expand with acceleration, 
i.e. when d 2 a/dt 2 = 0. Then one obtains 

^ 53.80 -I,o g (M^X). (12) 

Quantity (11) is related to the small slow roll parameters as follows [16] 

CM 

where 0^ and e refer to the moment of horizon exit and the end of inflation, correspondingly. 
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Introducing variable x — exp(^/2/3 0/Mp) we write down the slow roll parameters in R 2 - 
model (see e.g. [2]): 

Ml (V'\ 2 4 1 

^tA^W — ^' (15) 

(16) 

where prime denotes derivative with respect to scalaron field 0. Plugging (14), (15) into (13) 
and integrating then Eq. (13) one extracts slow roll parameters as functions of 

N=^N e + Xe -l, 



namely: 



4 1 /log(iV)\ 41 4 1 /log(iV)\ 4 1 ^f \og(N) \ 



3 iV 2 \ N 3 J ' ' 3 TV 3 7V 2 V iV3 / 3 ^ V Ar3 

Numerically, from the Friedman equation at the end of inflation, Xe ~ 4.6, and the relevant 
slow roll parameters are reasonably small, e e ~ 0.10, i] e ps 0.27, which justifies using of 
approximate relation (13) (see Ref. [16] for the exact relation). 

Tilts of the power spectra of scalar and tensor perturbations, (l-n s ) and tit, and tensor- 
to-scalar ratio r to the next-to-leading order in the slow roll parameters are given by [15, 16] 

2 9 >9 8 1 4.813 „ /log(iV)\ , x 

l-n s = 6 e -2r ? --r ? 2 + 0.374C 2 = -- + ^ + O(^^ r iJ , (18) 

Finally, substituting (9), (12) into (18)-(20) we obtain predictions for the cosmological pa- 
rameters. They are presented in Table 1, together with predictions for two other inflationary 
models exhibiting the same inflaton potential at inflationary stage: i? 2 -model with the Higgs 
field minimally coupled to gravity (1) and the so-called Higgs-inflation [8]. For these latter 
two models we refine the results of Ref. [7] for n s derived there to the leading order in slow 
roll parameters. The (absolute) error is expected to be of about 1/iV 2 ~ 10~ 3 . Indeed, our 
estimates of n s in these models deviate from those in Ref. [7] by this amount. 
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Model 


T rch , GeV 


n s 


r 


riT 


R 2 with conformal Higgs 


1.4 x 10 8 


0.9638 


0.0038 


-0.00047 


R 2 


3.1 x 10 9 [10] 


0.9644 


0.0036 


-0.00045 


Higgs-inflation 


6 x 10 13 [17] 


0.9664 


0.0032 


-0.00040 



Table 1: Next-to- leading order predictions for spectral parameters. 



For all the three models the interesting values of n s and r are well inside the region 
preferred by combined analysis of present cosmological data [18]. Nominally, the (absolute) 
error of the next-to-leading approximation \og(N)/N 3 = 10~ 5 is enough to ensure that it is 
possible to distinguish all three models by measuring the fourth digit in the values n s and r, 
parameter tit seems less promising. Provided the lower reheating temperature, predictions 
in our model deviate slightly stronger from those in the Higgs-inflation, as compared to the 
predictions in i? 2 -model with the Higgs field minimally coupled to gravity. Thus, future 
experiments, like CMBPol [19] with accuracy of 10 -3 in r and 0.0016 in n s , have better 
chance to distinguish the Higgs-inflation from our variant of R 2 , if any signal would point at 
the right ballpark. 

5 Gravity wave signals from inflation and from scalaron 
clumps 

In this model there are two sources of gravity waves: metric fluctuations at inflationary stage 
and scalaron clumps at late post- inflationary stage. Let us discuss them in order. 

The first source is inherent in any inflationary model. In the model under discussion 
it gives rise to tensor perturbations (gravity waves) with almost flat power spectrum after 
inflation: a deviation from flatness is characterized by index nx presented in Table 1. The 
total power in tensor perturbations is about 0.4% of that in scalar perturbations, see values 
of parameter r in Table 1. In the expanding Universe the perturbations, which length became 
smaller than horizon, start to evolve. Energy density of subhorizon tensor modes decreases 
with scale factor as radiation energy density, i.e. as 1/a 4 . Hence at post-inflationary stage 
when the Universe is dominated by oscillating scalaron with energy density scaling as 1/a 3 , 
the relative contribution of subhorizon gravity waves to total energy density drops as 1/a 
up to reheating, and later at radiation domination remains constant. Hence one expects a 
knee-like feature in the gravity wave spectrum to be observed at a frequency /* determined 
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by the horizon size at reheating i? re h- The latter is related to the reheating temperature 
through the Friedman equation, so that 

Hrch = 71 m p ■ (21) 



The relative contribution of tensor modes of conformal momenta 

k > /c rc h = 

<^reh 

are suppressed. At present /c re h corresponds to the physical momentum 



P* = — = — #reh • (22) 

a a 



By making use of entropy conservation we obtain 

To ( g* (T ) 



«o T rc h \g* (T, 



- reh y 



1/3 

(23) 



with effective degrees of freedom at present (T ) = 3.91 and g* (T reh ) = 106.75 [2]. Then 
substituting (23) and (21) into (22) we get for the present frequency, where the knee in 
gravity wave spectrum is expected, see Fig. 1, 

/„ = El. = 2.8 Hz ( ^ ^ . (24) 

J 2tt Vl-4xl0 8 GeVy V ; 

The second source of gravity waves is scalaron inhomogeneities: subhorizon modes at 
intermediate matter dominated stage grow proportionally to scale factor and have enough 
time to enter nonlinear regime before reheating [7] . Then, gravity waves can be produced at 
formation of scalaron clumps, at their subsequent merging and at final evaporation, see e.g. 
[21, 22]. The highest amplitude, hence the best chance to be observed, is expected from the 
latter process. Evaporation is the scalaron decays into relativistic SM particles. It is out-of 
equilibrium process yielding a non-vanishing transverse-traceless part of energy-momentum 
tensor feeding gravity waves. The typical frequency at production is about i/ rc h [21], which 
gets redshifted and presently coincides with /„ (24). The signal amplitude does not depend 
on the reheating temperature [21]. The estimate in [22] gives for a relative contribution of 
the gravity waves to the present total energy density Q gw ~4 x 10~ 13 £, where e < 1 is an 
efficiency factor which represents a measure of the asphericity of structure evaporation. The 
possible signal is shown in Fig. 1. 
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frequency [Hz] 



Figure 1: Energy density in gravity waves (in units of the present day critical density) Q gw as 
a function of frequency and the projected sensitivities of next-generation gravitational wave 
detectors: LIGO [20], BBO [5], DECIGO [6]. The picture shows the gravity wave signal 
from inflation (solid line) and from structure evaporation at reheating (red star). 

Similar signal is expected [7] in i? 2 -model with minimally coupled to gravity Higgs field. 
However, the typical frequency is by a factor T^/T re h higher. The same conclusion holds 
for position of the knee in gravity wave spectrum from inflation. One observes in Fig. 1 that 
those signals are either at the board of or out of reach of proposed experiments on searches 
for gravity waves. The signals in our model discussed above are right in the region available 
for investigation by future detectors like BBO [5] and DECIGO [6], see Fig. 1. This allows for 
independent test of our model: signatures in gravity waves pins down the value of reheating 
temperature. 
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6 Electroweak vacuum of the Higgs potential: lower 
limits on the Higgs boson mass 

Since we have modified the Higgs sector by introducing conformal coupling to gravity, the 
stability of the electroweak (EW) vacuum and whether the Universe ends up in it, given the 
cosmological evolution suggested in this paper, have to be investigated. The object under 
study is the Higgs effective potential, which in the unitary gauge "H T = (0, (h + v) /\/2) at 
large /i>u = 246.2 GeV reads 

V(h) = *&h'-±Rh*. (25) 

Here X(h) solves the renormalization group equation when h is replaced with renormalization 
energy scale [23, 24]. At large h selfcoupling X(h) may become negative, hence EW vacuum is 
metastable. Coefficient 1/12 in front of the second term in (25) is renorminvariant (neglecting 
graviton loops) [25]. For homogeneous, isotropic and flat Universe 

R = -12 H 2 -QH , (26) 

where dot denotes the time derivative. Thus, in the hot Universe: at radiation domination 
R = 0, at matter domination R = —3H 2 < 0, at present (dark energy domination), R < 0. 
Hence, the conformal coupling only enlarges the area of stability of the EW (i.e. widens the 
range of the allowed Higgs boson mass related to the selfcoupling constant as \/2Xv). 

Actually, metastability of the Higgs potential doesn't mean the theory is invalid. The 
weaker condition to require is the lifetime of the EW vacuum (with respect to tunneling 
and thermal decay) exceeds the Universe lifetime. For the usual case of minimally coupled 
Higgs the bound coming from tunneling is Mh > 111 GeV [26, 27] which is below the direct 
lower bound by ATLAS (115.5 GeV at 95% CL [28]). Since in our model the additional term 
RwH/6 stabilizes the potential, the condition from tunneling is fulfilled. 

The bound coming from the thermal decay depends on the reheating temperature [29, 
27]. The numerical calculation (3-loop /3-functions [14] and 0(a^), O (aa s ) corrections in 
matching of pole and running masses included [13]) yields a lower limit which is even below 
that from tunneling: 



M, 



therm 



M t - 173.2 GeV a, - 0.1184 

109.76 + — x 1.9 - x 0.14 

0.9 GeV 0.0007 



GeV . (27) 



Hence, given the direct limit from ATLAS, the problems with EW vacuum stability in our 
model are neither at hot stage of the Universe evolution, nor at present. 
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In our model the strongest lower limit on the Higgs boson mass comes from analysis of 
the Higgs field evolution from inflation to reheating. There at short periods of time At ~ 1/// 
the scalar curvature is positive because of oscillating scalaron (26), 

i?~-^(l-3cos(2^)) . (28) 

Then the Higgs potential (25) gets negative mass squared of order H 2 , which may stimulate 
escape to large values of fields and hence to wrong vacuum at late times. Let us analyze this 
process. 

Even if initially (classically) at origin, the Higgs field takes value h ~ H by the end of 
inflation, because of its quantum fluctuations, which get amplified when became superhori- 
zon. This estimate may be refined, e.g. by making use of the stochastic approach, see e.g. 
[30, 31]. The comoving probability P c (h,t) for the field to take value h at moment t obeys 
the Fokker-Planck equation 



dP c d 



H 3 dP c V'jh) 
8tt 2 dh ZH c 



(29) 



dt dh 

One can integrate it over h and introducing the Higgs correlators 

(h 2 ) = Jh 2 P c (h,t)dh , (hV'(h)) = JhV'(h)P c (h,t)dh 
cast it in the following form 

>> = £ - ■ < 30 > 

Neglecting the /i-dependence of A and treating perturbations as Gaussian quantities (A 
is reasonably small) we simplify (hV'(h)) ~ 3\(h 2 ) 2 + 2H 2 (h 2 ) and finally arrive at 

> = ^-|« 2 -f«. (3D 

Fluctuation (h 2 ) grows with time due to the first term and reaches maximum value when 
the r.h.s. of Eq. (31) is zero, that is 

^/WU^=^H. (32) 

After inflation the Higgs field starts from 2 (32) and evolves according to the classical 
equation of motion, 

h + 3Hh+ (-^R + \(h)h 2 j h = 0, (33) 



2 In fact the resulting lower limit on the Higgs boson mass (34) is (almost) insensitive to the numerical 
coefficient of order one in Eq.(32) in front of H . 
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where R(t) is given by Eq. (28). Given Eq. (28) the last term in parentheses in (33) is 
negligible at small t, and the field falls from initial value (32) as h oc t -1 / 3 oc 1/^/a. However 
at some moment the potential starts to dominate (A h 2 oc t~ 2 ^ 3 falls slower than R oc 2/9t 2 ) 
and if the corresponding value of h is above the maximum of the effective potential (25) 
then it rolls down to wrong minimum (towards large h). Numerical solution with X(h), 
evaluated by using 3- loop /3-functions [14] and O (a 3 ), O (aa s ) corrections in matching of 
pole and running masses [13], reveals that the Higgs field remains in small value region (and 
hence later evolves to the EW vacuum) provided the Higgs boson mass is above the following 
critical value 



M, 



crit 



M t~ 173.2 GeV _ a s - 0.1184 nn 

126.2 + — x 1.55 x 0.3 

0.9 GeV 0.0007 



GeV . (34) 



Uncertainties of estimate (34) are associated with errors in extraction of SM parameters 
M t and a s from experimental data (about 2 GeV at 65% CL) and unknown higher-order 
QCD-corrections to the effective potential (25) (about 1 GeV) [13]. Given these numbers, 
we conclude that obtained limit (34) does not contradict to the recent observation of the SM 
Higgs-like signals at LHC 

ATLAS [33] : M h = 126.0 ± 0.4(stat.) ± 0.4(syst.) GeV , (35) 
CMS [32] : M h = 125.3 ± 0.4(stat.) ± 0.5(syst.) GeV . (36) 

Therefore, being (minimally) extended with free scalar [34] or fermion [10] to serve as the 
dark matter, and with sterile neutrinos [10] to generate baryon asymmetry of the Universe 
and explain neutrino oscillations, the model we discussed becomes a phenomenologically 
complete, yet minimal, viable and testable model of particle physics. 
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